clear all; close all; clc;

R = 120; % Ом
C = 33e-12; % Ф
L = 100e-6; % Гн

omega0 = 1/sqrt(L*C); % Резонансная частота
f0 = omega0 / 2 / pi;
fprintf('f0 = %f MHz\n', f0/1e6);


T = 1 / (100*f0);
t = 0:T:(10 * 1/f0);
lt = length(t);

Uc = nan(1, lt);
dUc = nan(1, lt);
i = nan(1, lt);
di = nan(1, lt);

U0 = 10;
Uc(1) = U0; % Начальное напряжение конденсатора
i(1) = 0; % В начальный момент ток в цепи отсутствует
stdn = 3;
for k=2:lt;
    S(k) = randn(1,1) * stdn;

    Uc(k)=Uc(k-1)+i(k-1)*T/C;
 i(k)=(S(k)*T+L*i(k-1)-Uc(k)*T)/(L+R*T);
 Ur(k)=i(k)*R;
end;
figure(1);
plot(t,Ur, t,S);
figure(2);
plot(t, Ur)